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5BsinkgZ§y,  where  BQ=const.,  6B*const.,  and  X0*2ir/k0*const.  is  the  axial- wave¬ 
length  of  the  wiggler  field.  It  Is  assumed  that  the  beam  density  and  current 
ire  sufficiently  small  that  the  equilibrium  self  fields  can  be  neglected  in 
:omparison  with  B®.  In  this  context,  it  is  found  that  there  are  three  useful 
(and  exact)  invariants  (C_L,Cj1,Cz)  associated  with  single-particle  motion  in  the 
iquilibrium  field  Bg§z+5§.  These  invariants  are  used  to  construct  radially 
:onfined  Vlasov  equilibria  fj^Ci,.^,^)  for  an  intense  relativistic  electron 
aeam  propagating  primarily  in  the  z-direction.  Examples  of  both  solid  and 
lollow  beam  equilibria  are  considered,  and  it  is  shown  that  the  transverse 
higgler  field  can  have  a  large  modulational  influence  on  the  beam  envelope, 
iepending  on  the  size  of  SB/Bg. 
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ABSTRACT 

The  purpose  of  this  paper  is  to  develop  a  self-consistent  kinetic 
description  of  helically  distorted  relativistic  electron  beam  equilibria 
for  free  electron  laser  applications*  In  particular,  radially  confined 
equilibria  are  considered  for  a  helically  distorted  electron  beam  prop¬ 
agating  in  the  combined  transverse  wiggler  and  uniform  axial  guide 
fields  described  by  ^-BQjl^fi^-BQ^-fiBcoskQZ^-SBsinkQZ^y,  where  BQ* 
const*,  5B»const.,  and  XQ*2.7r/kQ- const,  is  the  axial  wavelength  of  the 
wiggler  field.  It  is  assumed  that  the  beam  density  and  current  are 
sufficiently  small  that  the  equilibrium  self  fields  can  be  neglected 
in  comparison  with  In  this  context,  it  is  found  that  there  are 

three  useful  (and  exact)  invariants  (Cx»C^,Cz)  associated  with  single¬ 
particle  motion  in  the  equilibrium  field  Bg£2+<$£*  These  invariants 
are  used  to  construct  radislly  confined  Vlasov  equilibria  f^CC^C^tC^) 
for  an  intense  relativistic  electron  beam  propagating  primarily 
in  the  z-dlrection*  Examples  of  both  solid  and  hollow  beam  equilibria 
are  considered,  and  it  is  shown  that  the  transverse  wiggler  field  can 
have  a  large  modulational  influence  on  the  beam  envelope,  depending 
on  the  size  of  6B/Bq. 
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I.  INTRODUCTION 

There  have  been  several  theoretical*  and  experimental^* ^ 
investigations  of  the  free  electron  laser  which  generates  coherent 
electromagnetic  radiation  using  an  intense  relativistic  electron  beam 
as  an  energy  source.  With  few  exceptions,  theoretical  studies  of  the 
free  electron  laser  instability  are  based  on  highly  simplified  models 
which  neglect  the  influence  of  finite  radial  geometry  and  beam 
kinetic  effects,  or  make  use  of  very  idealized  approximations  in 
analyzing  the  matrix  dispersion  equation.  The  purpose  of  this  paper 
is  to  develop  a  self-consistent  kinetic  description  of  helically 
distorted  relativistic  electron  beam  equilibria  for  free  electron  laser 
applications.  In  particular,  we  consider  radially  confined  Vlasov 
equilibria  for  a  helically  distorted  relativistic  electron  beam 
propagating  in  the  combined  transverse  wiggler  and  uniform  axial  guide 
fields  described  by  [Eqs.  (1)  and  (2)] 

«0-BoSr+{« 

■B0^“6Bcosk0Z^x'6B8ink0^y  ’ 

where  BQ*const. ,  6B*const. ,  and  XQ*2ir/kQ*const.  is  the  axial  wavelength 
of  the  helical  wiggler  field.  Both  solid  and  hollow  beam  equilibria 
are  considered,  and  it  is  assumed  that  the  beam  density  and  current 
are  sufficiently  small  that  equilibrium  self  fields  can  be  neglected 
in  comparison  with  B®.  Within  the  context  of  this  assumption  it  is 
found  that  there  are  three  useful  (and  exact)  invariants  (C^C^C  ) 
associated  with  single-particle  motion  in  the  equilibrium  field  Bq£z+6£< 
These  invariants  can  be  used  to  construct  radially  confined  Vlasov 
equilibria  f£  (CX,C ^»C2)  for  an  intense  relativistic  electron  beam 
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propagating  In  the  z-dlrection,  Including  the  Important  modulational 
Influence  of  the  transverse  wlggler  field  on  the  beam  envelope. 

The  three  exact  invariants  derived  in  Sec.  II  and  Appendix  A 
are  given  by  the  perpendicular  invariant  Cx  [Eq.  (4)], 

2  2  2eB0  2e 

c^"pr+pe+  <P2-*bmVb)+  d“  • 

the  helical  invariant  [Eq.  (5)], 

W  ^  <IVV,V+  • 

and  the  axial  invariant  C  defined  by  [Eq.  (6)] 


where  p  *  =  *s  t^ie  mec^ianical  momentum,  P0*r(p^-eBQr/2c) 

*  2 

is  the  canonical  angular  momentum,  Y^mc  *const.  is  the  characteristic 

directed  energy  of  the  electron  beam,  and  V^-const.  is  the  characteristic 

mean  axial  velocity.  In  the  expressions  for  Cx  and  C^,  we  have  subtracted 

the  constant  terms  proportional  to  Y^mV^  without  loss  of  generality. 

Moreover,  for  6B*0,  we  note  from  Eq.  (6)  that  the  axial  momentum  p  is 

z 

a  constant  of  the  motion,  and  that  Eqs.  (4)  and  (5)  yield  the  familiar 
2  2 

invariants,  p  +p  -const.  and  P  -const.,  for  a  charged  particle  moving 
r  “  o 

in  a  uniform  axial  field  Bq. 

The  striking  feature  of  the  present  analysis  is  the  fact  that  the 

exact  invariants  (C  ,C.  ,C  )  can  be  used  to  construct  helically 

i.  n  z 

distorted  relativistic  electron  beam  equilibria  f®  (C^.C^.C^)  °* 
experimental  interest  for  free  electron  laser  applications.  Generally 
speaking,  the  two  classes  of  relevant  beam  equilibria  can  be  charac¬ 
terized  as  (a)  solid  beam  equilibria,  and  (b)  hollow  beam  equilibria. 

For  example,  equilibrium  distribution  functions  of  the  form  [Eq.  (10)] 
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fb"Fl(C^"2Yb,nu)bCh)G(Cz)  ’ 

where  G(C  )  is  strongly  peaked  around  C  mV,  ,  correspond  to  a  solid 
electron  beam  with  non-zero  density  on  axis  (r«0).  On  the  other  hand, 
the  class  of  beam  equilibria  described  by  [Eq.  (11)] 

£®-F2(ei)«<ch-Vc<c,>  • 

2 

where  CQ=*-(eBQ/2c)RQ*const. ,  corresponds  to  a  slowly  rotating  annular 
electron  beam  with  characteristic  mean  radius  Rq.  Of  course,  it  is 
found  that  the  detailed  spatial  dependence  of  beam  equilibrium  prop¬ 
erties  (density  profile,  temperature  profile,  etc.)  depends  on  the 
specific  functional  form  of  G(CZ),  F2(Cx)  and  F^CCx-Zy^mw^C^) .  It 
is  also  found  that  the  radial  envelope  of  the  beam  can  be  strongly 
modulated  by  the  transverse  wiggler  field,  depending  on  the  size 
of  6B/Bq  [Sec.  III]. 

As  a  general  remark,  one  of  the  most  appealing  features  of  the 

3 

recent  free  electron  laser  stability  analyses  by  Davidson  and  Uhm, 
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Sprangle  and  Smith,  and  Bernstein  and  Hirshfield  for  a  relativistic 
electron  beam  with  uniform  density  and  infinite  cross-section  is  the 
fact  that  the  influence  of  the  transverse  wiggler  field  is  contained 
in  a  fully  self-consistent  manner  in  the  equilibrium  distribution 
function  f^-n^fi  (Px)6  (PyJGg^)  *  That  is,  when  carrying  out  the  stability 
analysis,  the  excited  electromagnetic  and  electrostatic  fields  are 
treated  as  small-amplitude  perturbations  about  a  self-consistent 
equilibrium  that  includes  the  full  nonlinear  influence  of  the  equilibrium 
wiggler  field.  We  believe  that  the  present  equilibrium  investigations 
form  an  important  first  step  in  formulating  a  self-consistent 
Vlasov  description  of  the  free  electron  laser  instability  that  Includes 
the  effects  of  finite  radial  geometry  and  also  correctly  incorporates 


the  nonlinear  influence  of  the  transverse  vlggler  field  on  the  bean 
equilibrium. 

The  organisation  of  this  paper  Is  the  following.  In  Sec.  II, 
we  discuss  the  basic  equilibrium  configuration  and  assumptions. 
Specific  examples  of  helically  modulated  relativistic  electron  beam 
equilibria  are  analyzed  in  Sec.  Ill,  both  for  solid  [Sec.  III. A]  and 
hollow  [Sec.  III.B]  electron  beams. 


II.  EQUILIBRIUM  CONFIGURATION  AND  BASIC  ASSUMPTIONS 


We  consider  the  class  of  helically  modulated  relativistic  electron 
beam  equilibria  propagating  in  an  externally  applied  magnetic  field 

«°-V*+s«  •  (1> 

where  Bq* const,  is  the  axial  magnetic  field,  and 

6B—  6Bcosk0zex-6Bsink0z£y  ,  (2) 

is  the  transverse  helical  wiggler  field  with  axial  wavelength  X^^ir/kg. 

In  the  present  analysis,  we  assume  that  6B«const.  ia  a  good  approximation 
over  the  radial  extent  of  the  beam.  In  cylindrical  polar  coordinates 
(r,0,z),  Eq.  (2)  can  also  be  expressed  as 

6B-6B  e  +6BaeA 
r^r  0^0 

— 6Bcos(e-k0z)dr+<5Bsin(e-k0z)£e  ,  (3) 

where  and  are  unit  vectors  in  the  r-  and  0-directions,  respectively. 

It  is  assumed  that  the  beam  density  and  current  are  sufficiently 

low  that  the  influence  of  the  equilibrium  self  electric  and  self 

magnetic  fields,  and  £g(/(0>  on  the  particle  trajectories 

0 

can  be  neglected  in  comparison  with  the  vxB  force  associated 
with  the  applied  magnetic  field  in  Eq.  (1). 

Within  the  context  of  the  above  assumptions,  there  are  three  useful  and 
exact  invariants  (Cj^C^C  )  associated  with  single-particle  motion  in 
the  equilibrium  field  These  invariants  can  be  used  to 

construct  cylindrical  Vlasov  equilibria  f£  (Cj^C^C^)  i°r  an  intense 
relativistic  electron  beam  propagating  in  the  z-direction,  including 
the  Important  module tional  influence  of  the  transverse  wiggler  field 


i 
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on  the  beam  envelope.  As  outlined  In  Appendix  A,  for  relativistic 
electron  motion  in  the  applied  magnetic  field  defined  in  Eqs.  (1)  and 
(2) ,  the  three  exact  invariants  are  given  by  the  transverse  invariant 
Cx  { Eq  #  (A.  12)  ] , 

2  2  2e 

C^Pr+pe+  (pz"YbmVb)+  d^  *  (4) 

the  helical  invariant  [Eq.  (A. 7)], 

W  ^  (pz-YbmVb)+  ^  rsin(e-k0z)  .  (5) 

and  the  axial  invariant  C  defined  by  [Eq.  (A. 11)], 

z 


V  ck 


f  eBoV  2e 


In  Eqs.  (4)-(6),  J>“(pr>Pg ,p2)“Ym^  is  the  mechanical  momentum,  P0« 

r(p  -eBnr/2c)  is  the  canonical  angular  momentum  associated  with  the 
0  U 

2  2  4  2  2  1/2 

axial  field  Bq,  ymc  «(m  c  +c  £  )  is  the  relativistic  electron  energy, 

2 

-e  is  the  electron  charge,  m  is  the  electron  rest  mass,  -const, 

is  the  characteristic  directed  energy  of  the  electron  beam,  and  V^* 
const,  is  the  characteristic  mean  axial  velocity.  Moreover, 
is  the  transverse  momentum,  and  can  be  expressed  as  £x*6£«- 

-pr6Bcos(0-kQz)+pQ<5Bsin(0-kQz) .  Note  in  Eqs.  (4)  and  (5)  that 
we  have  subtracted  the  constant  terms  proportional  to  Y^mV^  without 
loss  of  generality.  For  BQ-const.  and  6B-const.,  we  reiterate  that 
Cx,  C^,  and  Cz  are  exact  constants  of  the  motion  for  arbitrary 
wiggler  amplitude  6B. 

In  the  limit  of  zero  wiggler  amplitude,  6B-*0,  Eqs.  (4)-(6) 


reduce  to 


q22^  e®o 

Cx-Pr+Pe+  “d^  (pz-YbmVb)-const. , 


(7) 


(8) 


W  (Pz-Ybn>Vb)-con8t. . 

0 

C  “p  -const.  (9) 

z  z 

For  6B-0,  we  note  from  Eq.  (9)  that  the  axial  momentum  p  is  a 

z 

constant  of  the  motion,  and  Eqs.  (7)  and  (8)  yield  the  familiar 
2  2 

invariants,  p  +p  “const.  and  P  -const. ,  for  a  charged  particle 
moving  in  a  uniform  axial  field  Bq. 

There  are  two  classes  of  helical  beam  equilibria  f^CC^C^C^) 
of  experimental  interest  for  free  electron  laser  applications. 

Generally  speaking,  the  two  classes  can  be  characterized  as  (a)  solid 
beam  equilibria,  and  (b)  hollow  beam  equilibria.  For  example, 
equilibrium  distribution  functions  of  the  form 

fb"Fl(Cj-'2Ybm“bCh)G(C*)  »  (10) 

where  G(C  )  is  strongly  peaked  around  C  -yKmV,  ,  corresponds  to  a 

Z  Z  D  D 

solid  electron  beam  with  non-zero  density  on  axis  (r-0) .  In  Eq.  (10), 
ui^-const.  is  related  to  the  mean  angular  rotation  of  the  beam  [Sec.  Ill] 
On  the  other  hand,  the  class  of  beam  equilibria  described  by 

fb"F2(CA)6(Ch‘C0)G(Cz)  »  (ID 

2 

where  Cq— (eBQ/2c)RQ-const. ,  corresponds  to  a  slowly  rotating  annular 
electron  beam  with  characteristic  mean  radius  Rq.  Of  course,  the 
detailed  spatial  dependence  of  the  beam  equilibrium  properties  (e.g., 
density  profile,  temperature  profile,  etc.)  depends  on  the  specific 
functional  form  of  G(C^),  F^(CA-2y^m^C^)  and  F2(CX).  In  this  regard, 
specific  examples  of  helical  beam  equilibria  are  discussed  in  Sec.  III. 
As  a  general  remark,  it  is  found  that  the  radial  extent  of  the  beam 
can  be  strongly  modulated  by  the  transverse  wlggler  field,  depending 
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on  the  size  of  6B/Bq. 

While  it  is  true  that  beam  equilibrium  properties  can  be  calculated 
from  Eqs.  (10)  and  (11)  making  use  of  the  exact  invariants  in  Eqs.  (4)- 
(6) ,  there  are  useful  approximations  that  can  be  made  in  simplifying 
the  expression  for  C  [Eq.  (6)]  in  the  regimes  of  practical  interest 
for  free  electron  laser  applications.  We  now  discuss  these  approximations, 
which  will  be  used  in  the  detailed  equilibrium  examples  presented  in 
Sec.  III. 

First,  in  the  regimes  of  practical  interest,  the  characteristic 

transverse  momentum  j>A  of  a  beam  electron  is  small  in  comparison 

with  the  characteristic  directed  axial  momentum  Y.mV.'yp  ,  i.e., 

b  z 

IfcJ  <<  ^bmVb  •  02) 


To  assure  radial  confinement  of  the  beam  electrons,  the  present  analysis 
of  course  assumes  BqJ*0.  Making  use  of  the  inequality  |j>A|  << 

[Eq*  (12)],  the  striking  feature  of  the  present  analysis  is  that 
the  inequality  in  Eq.  (14)  is  easily  satisfied  even  for  large  wiggler 
amplitude  (6B%Bq,  say),  both  in  the  limits  of  weak  axial  field 
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2  2  2  2  2  2 
(o)c«4)q)  and  strong  axial  field  (wc>>Wq).  For  example,  if 

then  Eq.  (14)  reduces  to 


«c 

1»2  — 
“0 


Pi 

_ai _ 


*bmVb 


«§ 

Bn 


(16) 


which  is  readily  satisfied  in  the  parameter  regimes  of  experimental 

interest.  In  any  case,  within  the  context  of  the  inequalities  in 

Eqs.  (12)  and  (13),  the  exact  axial  invariant  C  defined  in  Eq.  (6) 

z 

can  be  approximated  by 


C  *p  - 
z 


ck0'Ei‘,SE 

3> 

Pz  ckft 


(17) 


where  we  have  taken  the  positive  square  root  in  Eq.  (6),  which 

corresponds  to  C  and  p  having  the  same  polarity.  Equation  (17) 
z  z 

is  an  excellent  approximation  to  C  in  the  parameter  regimes  of  practical 

z 

interest  for  free  electron  laser  applications.  Moreover,  for  G(C  ) 

z 

strongly  peaked  about  with  characteristic  half-width 

ACz<<YbmVb>  t*le  denominator  on  the  right-hand  side  of  Eq.  (17) 
can  be  approximated  by  pz-eBQ/ckQ=Y^mV^-eBQ/ckQ.  Equation  (17) 
then  reduces  to 


wc 

^z  ^z  w  -w  'VL  Bn  9 
0  c  0 

which  is  the  approximate  form  of  the  axial  invariant  used  in  Sec. 


(18) 

III. 


III.  EXAMPLES  OF  HELICALLY  MODULATED  RELATIVISTIC 


ELECTRON  BEAM  EQUILIBRIA 


There  is  clearly  a  wide  variety  of  helically  modulated  relativistic 
electron  beam  equilibria  that  can  be  analyzed  within  the  context  of  the 
equilibrium  formalism  described  in  Sec*  II*  The  relevant  choice  of 
(Cj^C^C^)  course  depends  in  detail  on  injection  geometry, 
beam  quality,  etc.  For  our  purposes  here,  we  consider  two  simple 
examples  of  solid  and  hollow  beam  equilibria  that  clearly  illustrate 
the  strong  influence  of  the  transverse  wiggler  field  in  modulating 
the  beam  envelope.  In  both  cases,  for  simplicity,  we  assume  that 
the  axial  distribution  G(Cz)  is  cold,  i.e., 


G(Ci).«(Crybmvi) 


to 

„  c  ^ 

p  -y,  mV, - p,  •  — 

Yz  b  b  u)A-a)  Ba 

U  c  u 


(19) 


where  use  has  been  made  of  Eq.  (18).  The  analysis  can  be  extended  in  a 
straightforward  manner  to  distribution  functions  G(CZ)  with  a  small 
spread  about  C  -y.mV,  . 

Z  D  O 


A.  Solid  Beam  Equilibria 


As  an  example  of  a  helically  modulated  solid  electron  beam,  we 
consider  the  case  where  f^(CA,C^,Cz)  has  the  form  [Eq.  (10)] 


fb“  ^  «<CA-*Tb«sV2^) 


(20) 


x6(WV  * 


where  and  are  defined  in  Eqs.  (4)  and  (5),  Cz  is  defined  in 
Eq.  (18),  and  ng,  and  TA  are  positive  constants.  Equation  (20) 
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is  a  straightforward  generalization  of  the  uniform  density  beam 
0  2  2 

equilibrium  fb-(no/ff)6(pr+p0-2ybma>bPe-2YbmTx)6(pz-YbinVb)  previously 

g 

discussed  by  Davidson  for  the  case  6B*0.  Evaluating  Cx-2Y^mo)t)Cji 
for  and  rearranging  terms ,  we  find  that 


V“b 


Pr+  ckQ  U0-o)c 


/OIq-U. 

6B  +  p.-y,  mw.r+  I - 

r  6  b  b  ck0  I  O)0-O)c 


+Ybm2\|;(r,e-k0z)  , 


where  w^eBQ/y^mc,  61^— 6Bcos  (S-k^z) ,  6Be*6Bsin(0-kQz) ,  and 

the  effective  potential  ^(r,e-kQz)  is  defined  by 

2  2 

*(r,  0-kQz)«  )  r 

/  \  /  \2  2  0 

..  /V“t>V  6B  0-“bY“c/««\V 

We  now  evaluate  various  macroscopic  properties  of  physical 
interest  for  the  choice  of  distribution  function  in  Eq.  (20).  For 
example,  making  use  of  Eqs.  (20)  and  (21),  it  is  readily  shown  that  the 
equilibrium  density  profile  n°(r,e,z)mJd3p  f^(CltC^,Cz)  corresponds 
to  a  constant-density  beam  with  sharp  radial  boundaries,  i.e.. 


nb(r,e-kQz)- 


nQ  ,  t(i(r,e-k0z)<2Tx/Ybm  , 

0  ,  lp(r,e-k0z)>2fx/ybm  , 


where  \p  is  defined  in  Eq.  (22).  Moreover,  making  use  of  |£x|«ybmVb, 
the  mean  equilibrium  velocity  of  the  electron  beam  can  be  approximated  by 
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Vb(^)“t/d3p(^/Ym)fb1/(^d3p  fb) 


(24) 


“  7^  ^  «  *!>'</*  fS>  • 

b 

and  the  transverse  temperature  profile  can  be  approximated  by 

T2b<JE>-  T  <{'*3P[<pr-<pI.>Hvr-<vr>)+<pe-<pe>)<ve-<.9>]fJ)/(/d3p£j) 


-2V  (/»3P I  cpr- <Pr»2+(pe-  <pe»2 If °)  / (/d3pfS>  , 

where  <pj>°(|d^ppjf^)/ (Jd\f®) ,  and  use  has  been  made  of  y-Y^ 
Substituting  Eqs.  (20)  and  (21)  into  Eq.  (24),  the  mean  macroscopic 
velocity  components  of  the  electron  beam  can  be  expressed  as 

Vjk(r,.-V)-Vb  co»(9-V>  • 


(25) 


(26) 


V"  (r,B-k  z)=ol  r-V  -£(  -2-^ )  |£  sin(0-knz) 
0b  0  b  b  wn  \  Wq-u>c  /  Bq  0 


and 


Vzb(r*e-k0z)=Vb  ’ 


(27) 


(28) 


in  the  region  of  configuration  space  where  ni;(r,e-knz)  is  non-zero. 

b  u 

In  obtaining  Eqa.  (26)  and  (27),  use  has  been  made  of  6Br*-6Bcos(e-kQz) 
and  6B^»6Bsin(6-kQz) ,  and  the  contribution  to  of  order 

has  been  neglected  in  obtaining  Eq.  (28).  Finally,  making  use  of  Eqs. 
(20),  (21),  and  (25),  the  transverse  temperature  profile  can  be 
expressed  as 

i|>(r,e-k0*)^ 

T 

'0 


-  t  i|>(r,e-kQz)i 
*2b  (r ,  6-k0z)«Tx  1 - ^ - 


(29) 


in  the  region  of  configuration  space  where  n^(r,6-knz)  is  non-zero, 

b  0 

i.e.,  where  i|»  <  2Tx/y.®  =  [Eq.  (23)]. 
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Evidently,  from  Eq.  (23),  the  outer  radius  R^(6-kgZ)  of  the  electron 

beam  is  determined  self-consistently  from  the  solution  to  the  "envelope" 

2 

equation  iKR^,0-kQz)=2TA/Y^m  =  Vq.  Making  use  of  Eq.  (22),  we  find 
that  the  physically  allowed  value  of  R^B-k^z)  is  given  by  the  expression 


o  2  2  1/2 

Rb(e-k0z)=4Rbsin(e-k0z)+[(ARb),dsin  (e-k0z)+Rj]‘1' 


(30) 


where  R_  and  AR,  are  defined  by 

0  D  \2  2 


*o- 


and 


2j^%\  !c/6B\V 
0  yvBcj 

KVV 


kQAR  =  — - — 

0  'b  u)Q-toc  Bq 


(31) 


(32) 


Moreover,  the  equilibrium  density  profile  in  Eq.  (23)  can  be  expressed  as 

n0  *  0  <  r  <  V0-koz)  » 

0  ,  r  >  Rb(e-kQz)  . 


n^(r,e-koz)a 


(33) 


In  addition,  after  some  straightforward  algebraic  manipulation  that 
makes  use  of  Eqs.  (29)-(32),  the  transverse  temperature  profile  can 
be  expressed  as 


0  (1-  R.  )  (1+  R.  "  Rfa  8in(0_ko*) 

T?b(r,e-k0z)=T”b  - ^  - -  (34) 

|l - jj-=-  sin(e-k(jz)J 

where  R^e-kgZ)  is  defined  in  Eq.  (30),  and  T^b  is  the  maximum,  on-axis 

(r>0) ,  transverse  temperature  defined  by 

2  2 


^2f<3£)3(crc- 


(35) 


As  a  simple  reference  equilibrium,  for  6B*0,  we  note  from 
Eqs.  (30)  and  (31)  that  the  beam  radius  R^  reduces  to  the  familiar 
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8  2  1/2 

result  »  an<*  *kat  AR^*Ot  which  corresponds  to 

zero  modulation  of  the  beam  envelope.  In  addition,  for  6B»0,  it  follows 
that  vJb-0  [Eq.  (26)],  vj^r  [Eq.  (27)],  and  T^-T^l-r^R2) 

[Eqs.  (34)  and  (35)].  That  Is,  for  zero  viggler  amplitude,  Eqs.  (30)- 

g 

(35)  reduce  to  the  familiar  results  corresponding  to  an  electron  beam 

with  uniform  density,  constant  beam  radius,  rigid-rotor  angular  velocity 

profile,  and  parabolic  transverse  temperature  profile. 

In  the  case  where  6B^0,  we  note  from  Eqs.  (30)-(32)  that  the  outer 

envelope  of  the  beam  is  modulated  by  the  transverse  wiggler  field. 

2 

Moreover,  for  existence  of  the  equilibrium,  >  0  is  required 
[Eq.  (31)]  so  that  must  lie  in  the  range 

0  <  %  <  wc  »  (36) 

for  radially  confined  equilibria.  From  Eq.  (30),  the  maximum  (R^) 
and  minimum  (R^)  radial  excursion  of  the  beam  envelope  occurs  for 
6-kQZ-(2n+l)7r/2,  n=0,  ±1,  ±2,...,  with  R^  and  defined  by 

and  (37) 

1  I  **-t  Z+R§  ] r/ 2  , 

where  Rq  and  AR^  are  defined  In  Eqs.  (31)  and  (32),  respectively. 
Defining  the  average  beam  radius  by  Rb“(R^-R^)/2  gives 

Rh-KAR^+R2]172  .  (38) 

Moreover,  defining  the  peak  amplitude  of  the  radial  modulation  by 
AR^O^-R^)/!  gives 
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AVlARbl 


(39) 


1_ 

0) 

c 

6B 

k0 

o 

3 

1 

o 

3 

B0 

We  now  parameterize  the  equilibrium  properties  within  the  context 
of  the  basic  assumptions  enumerated  in  Sec.  II.  First,  the  transverse 
(r,6)  motion  is  assumed  to  be  nonrelativistic  [Eq.  (12)].  Making  use  of 
| V^b I »  lV0bl  <<  c*  we  ^rom  E^8,  (26)  and  (27)  that  the  quantities 

o>c»  6B,  etc.  are  required  to  satisfy 

|Vb  “c  W  6b| 


c  “0  (u,0'“c)  B0 


«  1  , 


(40) 


and 


2-2 


«  1  . 


(41) 


where  uA  and  u  are  far  removed  from  cyclotron  resonance  [Eq.  (14)]. 
u  c 


For  a  weak  axial  field  with  o>q  >>  a>c  >  ,  Eq.  (40)  reduces  to 

«  1  ,  (42) 


Vb  “c  £B 


c  u0  B0 


which  is  easily  satisfied  even  for  moderate  wlggler  amplitudes  with 
|  6B/Bq  |  \  1.  On  the  other  hand,  for  a  strong  axial  field  with  Uq  «  &>c, 
Eq.  (40)  reduces  to 

I V. 


_b  6B 

c  BA 


«  1. 


(43) 


which  requires  a  small  wiggler  amplitude  with  |6B/Bq|  <<  1.  Here  we 
have  assumed  >  o^.  Finally,  for  present  purposes,  we  estimate 
o£r£/c2  w^R^/c2  £  “^vq/ WjjUg c^ [ Eq .  (31)].  The  Inequality  in  Eq.  (41) 

then  becomes 


“b  V0 


u  2 
c  c 


«  1  . 


(44) 
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which  is  easily  satisfied  since  (u^  <  u>c  [Eq.  (36)]  and  the 

2  2 

transverse  thermal  motion  is  nonrelatlvistic  (Vq  «  c  ). 

It  is  also  interesting  to  determine  the  characteristic  size  of 

the  normalized  maximum  radial  excursion  of  the  beam  envelope. 

From  Eq.  (39),  for  <i)Q  >>  o»c,  we  find  kg/TR^  =  |  (6B/Bq)|  «  1 

[Eq.  (42)]  for  weak  axial  magnetic  field.  Moreover,  for  wA  «  u  , 

u  c 

we  find  k^AR^  -  |6B/Bq|  <<  1  [Eq.  (43)]  for  strong  axial  magnetic 
field.  That  is. 


‘A 


«  1 


(45) 


follows  consistently  from  the  assumption  of  nonrelatlvistic  transverse 
motion  in  both  the  weak  field  [Eq.  (42)]  and  strong  field  [Eq.  (43)] 
regimes . 

We  also  determine  the  characteristic  size  of  k^R^.  Estimating 

R0  %  ^“^“c  ^Eq*  8iv*s 

2  2 

2_2  “0  ^ 
oV  “b“c  V*  ' 

D 

where  uo“R0^b  an<*  “b  <  “c  The  regime  of  most  interest 

2  2 

experimentally  is  kgR^  >  1,  which  requires  sufficiently  slow  beam 
rotation  satisfying  [Eq.  (46)] 


(46) 


<*>■ 


2  2 
b  u0  v0 


u> 


2  y2  * 

u  vr 

c  b 


(47) 


As  a  numerical  example,  in  Fig.  1  we  show  a  plot  of  the  normalized 
radius  k^R^  of  the  beam  envelope  [Eq.  (30)]  versus  S-k^z  for  several 
values  of  6B/Bq  and  the  choice  of  equilibrium  parameters  “o"^“c* 
u»k*0.1  u>c»  Vq/c3“0.04,  and  V3/c3«8/9.  In  this  case, 
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W  i  ^  >  “•> 

and 

1c2rJ-8[1+|  (19.5/9)2(«B/B0)2]  (49) 

follow  from  Eqs.  (31)  and  (32).  Substituting  Eqs.  (48)  and  (49)  into 
the  expression  for  R^  in  Eq.  (30)  gives  the  results  shown  in  Fig.  1 
for  6B/BQ-l/4  [Fig.  1(a)]*  6B/B0«0.5  [Fig.  1(b)]  and  6B/BQ-1  [Fig.  1(c)]. 
Ue  note  from  Eqs.  (38),  (39),  (48),  (49),  and  Fig.  1  that  the  normalized 

V*b 

increase  from  (k0Rb,k0ARb)*(2.83,0)  for  dB/BQ-0,to  (kQ^fk^AR^)- 
(5.185,  0.333)  for  6B/Bq=1. 

We  conclude  this  section  by  emphasizing  that  the  theoretical  model 

developed  here  can  be  used  to  calculate  the  equilibrium  properties 

of  a  helically  modulated  solid  electron  beam  propagating  in  an 

equilibrium  magnetic  field  prescribed  by  Eqs.  (1)  and  (2)  for  a  broad 

2  2 

range  of  system  parameters  u)q,  vq/c  »  6B/BQ,  etc.  The  major 

approximations  relate  to  the  assumptions  that  the  transverse  particle 
motion  is  nonrelativistic  [Eq.  (12)]  and  that  the  axial  motion  is 
far  removed  from  cyclotron  resonance  [Eq.  (14)]. 


average  beam  radius  k^R^  end  maximum  modulation  amplitude 


B.  Annular  Electron  Beam 


As  an  example  of  a  helically  modulated  annular  electron  beam, 
we  consider  the  case  where  ^vCA,C^,C2)  has  the  form  [Eq.  (11)] 

n  nn^n^Kmvn  o  ?  o  2  o 

- “  U[(Cx-Y‘m  vj)  -c2]6(Ch.C0)6(Cz-YbmVb)  , 


where  Cq--cBqRq /2c® const . ,  and  n^,  Rq,  v^,  and  e  are  positive 

constants.  In  Eq.  (50) t  U(x)  is  the  Heaviside  step  function  defined 

by  U(x)*l  for  x  >  0,  and  U(x)®0  for  x  <  0.  Evaluating  C^-Cq  and 
2 

Cx-(Yb®VQ)  at  cz*Y^mV^»  and  rearranging  terms,  we  find 

‘Wc.-vvh  k  “"<‘i-v>)n 


where 


kT^rlr'M<e-V)'';-F(t-9-V)  • 

0  0  c  0 


ICl*  <Yb»V  2  ]  Cz.Yb„Vb-P;2^;2- (I'b"*  2  W02 


viSX  bA(?J  • 


U) 

p’-p  -y.  mV.  — - — —  cos(8-knz)  , 
r  r  b  b  <*>q~<*>c  bq  0 


w  fiB 

Pfl“Pft+YK®VK - —  8in(e-kftz)  , 

8  0  b  b  u>q— u>c  Bq  0 


and  F(r,8-kQz)  is  defined  by 


F«y.mV.  fi  —  kn(r^-R?)+r  —  — — —  ~  sin(8-kftz) 
Tb  b  2  uu  0  0'  hu-m  Brt  O' 


i] 


It  is  convenient  to  introduce  the  quantities  f  and  g  defined  by 


~~  ~  I5-  sin(e-kQz)  , 
k0  “0  “c  B0  0 


g°  r~  -■  c(.  cos(e-knz)  . 
k0  “o  “c  B0  0 

2 

Furthermore,  evaluating  Cx-(ybmv0)  at  Ch*CQ  and  c2*YbmVb,  we  find 
from  Eqs.  (51),  (52),  and  (55)  that 


tCx-(Yb-*0)  1C 

Z  D  D 


VC0 


-((V>vo)2" 


(r+f)2+g2J 


where  is  defined  by 

0  „  ..  .  “c  «b  f(w-f)V-F/rbmv  k  * 

Pr-Tb„Vbco.(8-V)— -(  (r+f)2+g2 - j  ■  «, 

As  a  simple  example,  we  consider  the  limit  of  zero  perpendicular 
energy  spread  in  Eq.  (50).  Making  use  of  the  identity 

«(Cx-Yb®2vJ)“lim^- U[(CA-yJm2VQ)2-e2]  ,  ( 5< 

we  find  from  Eqs.  (50) ,  (51),  and  (56)  that  the  electron  density 
profile  n^(r,e,z)*/^a>dpr/^codp0/^oi(dp2  fb^Ci,Ch,Cz^  can  be  exPressed  as 


nb(r,e-k0z)* 


n0Vbmv0  r 

[(r+f)2+g2jl/2  J J 


j"  dp^(p"2-i».)  , 


where  p£-{l+[g/(r+f)]2)2^2(pr-p®) ,  and  is  defined  by 


r,e-k0z).( 


,  •  v2  F“ 

(Ykmvl) - 5 — 5-  . 

b  °  (r+f)2+g2J 


Carrying  out  the  p£  integration  in  Eq.  (59),  the  density  profile 
can  be  expressed  as 


noVo 
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n?(r,6-kftz>-  {v'2[(r+f)V]><P/Yb»)2>1/2 


♦  >  0  , 


•b'*’w  -0“'-  \  u  b  '  '  (61) 

0,  otherwise  . 

In  the  subsequent  analysis  It  Is  useful  to  Introduce  the  effective 
radial  variable  R  defined  by 

2 


8i"<9-lo°) 


_ 6B 

k2  l  V“c/  '"0 


(ir)  U-sin2(e-k0z)]  , 


and  the  constant  quantities  r^,R,R^,  and  defined  by 


(62) 


CL  fa)  * 

c 


(63) 


R2- 


c  0  c  ® 


K-h+<A4)U2 


(64) 

(65) 


V-'l+(52+il)1/2 


(66) 


Making  use  of  Eqs.  (61)-(66),  and  the  definition  of  i»i(r,e-k0z) 
in  Eq.  (60),  we  find  that  the  density  profile  in  Eq.  (61)  can  be 
expressed  in  the  equivalent  form 


n^r.e-kQ*)- 


n0R0 

[(r^2-r2)(r2-r^2)]1/2 


0  ,  otherwise  . 


.<  *  .<  «b 


(67) 


From  Eqs.  (61)  and  (67),  note  that  the  normalization  of  n^  is  such  that 
nb“nQ-const.  for  r-RQ  and  6B-0.  In  addition,  note  from  Eq.  (67)  that 
the  beam  density  is  singular  at  the  inner  (R-R^)  end  outer  (R»r£) 
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boundaries  of  the  electron  beam.  This  is  a  consequence  of  the  singular 

form  for  the  perpendicular  energy  distribution  function  assumed  in 

Eq.  (58).  For  small  but  finite  e,  the  beam  density  remains  finite 

o 

and  varies  smoothly  throughout  the  annulus.  We  emphasize,  however. 


that  the  beam  density  profile  in  Eq.  (67)  is  integrable  in  the  sense 

r  2tt 

i 

that  the  number  of  electrons  per  unit  axial  length 


-  v/  -  c 


dr  r  % 


finite.  In  any  case,  the  equilibrium  example  in  Eq.  (67)  is  adequate 
for  present  purposes  of  describing  the  helical  modulation  of  the  inner 
and  outer  beam  envelopes  by  the  transverse  wiggler  field. 

Making  use  of  Eqs.  (62),  (65),  (66),  and  (67),  the  outer  boundary 
(r*r*)  and  inner  boundary  (r=r£)  of  the  annular  electron  beam  are 
determined  from 


r*»-  - - - —  -jp  sin(0-knz) 

b  k0  “o" “c  B0  0 


r2^_2.1/2,2  1 


[rT+(R  +rr')  /  ]  -  ^ 


0  c 


llj)  [1-sin  (0-koz)1 


(68) 

1/2 


and 


r.  ■-  - - —  sin(0-knz) 

b  k0  “0'“c  B0  0 


52j._2»1/212  1 


[~rT+(R  +rr')  '  ]  -  ^ 


sin  (0-kQz)] 


1/2 


(69) 


where  r^v^/h^  m(^  B  *s  defined  In  E9*  (64)*  In  the  special  limiting 

case  where  6B*0,  it  follows  from  Eqs.  (63)  and  (64)  that  rL*vo^“c 

_  q 

and  R*Rq,  and  Eqs.  (68)  and  (69)  reduce  to  the  familiar  results 


and 


/  /1B2.  2v  1/2 

(rb)6B-0”(R0+rL)  ~rL  * 


(70) 


is 


(71) 


A*- 
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That  is,  for  6Bm0,  the  average  beam  radius  and  beam  thickness  are 
constant,  with  R^= (r^+rb)/2=(RQ+rL>  '  =  Rq  (for  r£  «  R*) ,  and 

iRb®VV2rL' 

For  the  choice  of  equilibrium  distribution  function  in  Eq.  (50) 
with  t  0+,  Eqs.  (68)  and  (69)  give  very  precise  predictions  for  the 
inner  and  outer  boundaries  of  the  annular  electron  beam  for  6B^0, 
within  the  context  of  the  assumptions  enumerated  in  Sec.  II. 

In  the  remainder  of  this  section,  we  make  use  of  Eqs.  (68)  and  (69) 
to  examine  the  helical  modulation  of  the  beam  envelope  in  the  general 
case  where  6 B  is  non-zero,  assuming  in  addition  that  the  dimensionless 
quantities 


r£/R2  <  1,  (72) 

and 


can  be  treated  as  small  parameters.  In  Eq.  (72),  R  is  defined  in 

Eq.  (64),  and  r^=2vQ/u)c,  w^ere  vq  is  defined  In  Eq.  (53),  in  the 

2  2 

regime  of  experimental  interest  RqRq  >  !•  In  this  regard,  we  note  that 

the  inequality  in  Eq.  (73)  is  readily  satisfied  even  for  moderate 

higgler  amplitudes  with  ) 6B/Bq |  %  1.  This  is  true  both  for  weak 

axial  fields  (co^  «  oJq)  and  strong  axial  fields  (a>c  >>  w^).  Taylor 

expanding  the  expressions  for  rg  in  Eqs.  (68)  and  (69),  and  retaining 
2  —2  2 

terms  to  order  (r^/R  )  and  (<$B/Bq)  ,  we  obtain  the  approximate 
expressions 


+ 


R  +  R 


1  Mc  6B 
kQR  “<T“c  B0 


sin(e-kgZ) 


(74) 


—  isa 

*ll  iTlTdfitei 


Evidently,  from  Eq.  (74),  the  beam  thickness  AR^r^-r^  is  given  by 


-  2  2 

correct  to  order  (rL/R)  and  (6B/Bq)  .  In  Eq.  (75),  v^  is  defined 
in  Eq.  (53)*  Moreover,  the  average  beam  radius  R^«(r^+r”)/2  is  given 


5  1- 


-  —  sin(6-knz) 

.  r  wn-u>  Bn  0 

kgR  0  c  0 


+  *  '  !  '  t2j2  (v«c)(»o)  [1'el'‘2<6‘V)1|  ’ 

-  2  2 

correct  to  order  (r^/R)  and  (6B/Bq)  .  We  note  from  Eq.  (76) 
that  the  leading-order  modulation  of  the  average  beam  radius 
is  linear  in  6B.  Therefore,  making  use  of  Eq.  (64),  the  expression 
for  R^  in  Eq.  (76)  can  be  approximated  by 

R»R  ]  1 - —  — ~  sin(0-kQz)  |  ,  (77) 

15  (  k0R  “0  “c  B0  °  1 

correct  to  lowest  order  in  6B/Bq. 

As  a  numerical  example,  in  Fig.  2  we  show  a  plot  of  the  normalized 
average  beam  radius  k^R^  (dashed  curve)  versus  G-k^z  obtained  from 
Eq.  (77)  for  several  values  of  6B/Bq  and  the  choice  of  equilibrium 
parameters  v2/c2*0.04,  kQR0*3/2,  k0rL“k0v0/a)c«l/20f  yb*3  and 


Vb/c2*8/9.  In  this  case 


k0*b"k0*'  1 . 


follow  from  Eqs.  (53),  (75),  and  (77).  Here,  kQR«kQRQ[l+(28/81) (6B/Bq)  ] 
follows  from  Eq.  (64).  Equations  (78)  and  (79)  yield  the  results  shown 


2,1/2 
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in  Fig.  2  for  6B/BQ«l/4  [Fig.  2(a)],  6B/BQ-0.5  [Fig.  2(b)]  and  6B/BQ-1 

[Fig.  2(c)].  From  Eq.  (79),  we  define  the  peak  amplitude  of  the  radial 

modulation  by  k^AR^  ■  | (1/3) (6B/Bq) | .  We  note  from  Fig.  2  and 

Eqs.  (78)  and  (79)  that  the  normalized  beam  thickness  k^ARfe  and 

maximum  modulation  amplitude  kgAR^  increase  from  (kgAR^,  kgAR^)* 

(0.10,  0)  for  6B=0,  to  (k^R^  kQARb)-(0.19,  0.33)  for  6B/Bq-1. 

In  concluding  this  section,  we  emphasize  that  the  theoretical 

model  developed  here  can  be  used  to  calculate  the  equilibrium  properties 

of  a  helically  modulated  annular  electron  beam  propagating  in  the 

equilibrium  magnetic  field  BQ^-fiBcoskQZ^-fiBsink^zjl^  for  a  broad 

2  2 

range  of  system  parameters  a>c,  u>q,  Vq/c  ,  6B/Bq,  k^R^,  etc.  As  in 
Sec.  II. A,  the  major  approximations  relate  to  the  assumptions  that 
the  transverse  electron  motion  is  nonrelativistic  [Eq.  (12)]  and  that 
the  axial  motion  is  far  removed  from  cyclotron  resonance  [Eq.  (14)]. 


IV.  CONCLUSIONS 


In  this  paper  we  have  developed  a  self-consistent  Vlasov  descrip¬ 
tion  of  helically  distorted  relativistic  electron  beam  equilibria 
for  free  electron  laser  applications*  In  particular  we  have  considered 
radially  confined  equilibria  for  a  helically  distorted  relativistic 
electron  beam  propagating  in  the  combined  transverse  wiggler  and  uniform 
axial  guide  fields  described  by  Eqs.  (1)  and  (2).  Assuming  that  the 
beam  density  and  current  are  sufficiently  small  that  equilibrium  self 
fields  can  be  neglected  ip  comparison  with  it  is  found  that 

there  are  three  useful  and  exact  invariants  (C^C^tC^)  associated  with 
single-particle  motion  in  the  equilibrium  field  configuration.  These 
invariants  are  used  to  construct  radially  confined  Vlasov  equilibria 
f£(Cx,VCx)  for  an  intense  relativistic  electron  beam  propagating 
primarily  in  the  z-direction.  Specific  examples  of  solid  [Eq.  (20)] 
and  hollow  [E ; .  (50)]  beam  equilibria  are  analyzed  in  Sec.  Ill,  and  it 
is  shown  that  the  transverse  wiggler  field  can  have  a  large  modulational 
influence  on  the  beam  envelope,  depending  on  the  size  of  6B/Bq.  We 
believe  that  the  present  equilibrium  investigations  form  an  important 
first  step  in  formulating  a  self-consistent  Vlasov  description  of 
the  free  electron  laser  instability  that  both  includes  the 
effects  of  finite  radial  geometry  and  correctly  incorporates  the  non¬ 
linear  influence  of  the  transverse  wiggler  field  on  the  beam  equilibrium. 
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FIGURE  CAPTIONS 

Fig.  1  Plot  of  [Eq.  (30)1  versus  0-kg*  for  <i»Q*4tttc,  (^-0.1  uc> 

Vq/c2“0.04,  Yb-3,  V2/c2-8/9  and  (a)  6B/B0-l/4,  (b)  «B/B0-l/2, 
and  (c)  6B/Bq*1. 

-  2  2 
Plot  of  kpR^  [Eq.  (77)1  versus  e-kQz  for  u0-4«c,  vQ/c  -0.04, 

yb-3,  V^/c2-8/9,  k0R0-3/2,  k0rJ-l/20,  and  (a)  6B/BQ-l/4, 

(b)  6B/Bq-1/2,  and  (c)  6B/BQ-1. 


Fig.  2 
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APPENDIX  A 

SINGLE-PARTICLE  CONSTANTS  OF  THE  MOTION 

In  this  appendix,  we  consider  the  motion  of  an  individual  electron 
in  the  applied  magnetic  field 

fc0“BoSz'6Bco8koz£x'6B8lnk°zSy 


mBo«z‘6Bcos(e'V)^r+6B8in(e'koz)^e  * 

where  B^- const.  is  the  uniform  axial  field,  6Jg“-6Bcos(6-kQz)j£r+6Bsin(6-k.Qz)j|e 

is  the  transverse  helical  wiggler  field  expressed  in  cylindrical 

polar  coordinates,  6B-const.  is  the  wiggler  amplitude,  and  X^-Eit/kg  is 

2  22  1/2 

the  wavelength.  Making  use  of  y  '-(1+jj'  /me)  "const.,  the  equations 
of  motion  for  the  particle  orbits  r'(t'),  e'(t’)  and  z'(t')  can  be 
expressed  as 


t'.(r'  .’♦*<&-)  (#)  )  ■  &  ♦  *  cos(9'-k0i')  &  (*.2) 

t'«  ^7f  -  -  “•<9'-V')c'  5T-  •  <A-3) 


Conservation  of  energy  (y' -const.)  of  course  implies  that 


P^2+Pg2+P*2“const. , 


(A.  4) 


where  p^-y’mdr'/dt' ,  p^-y'mr'deVdt'  and  p^-y'mdz'/dt' .  Moreover, 
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multiplying  Eq.  (A. 2)  by  r*  and  making  use  of  Eq.  (A. 3)  gives 
eB, 


(r’(pe-  i~r r')j-  r’co8<e,-ko*,)^- 


(A.  5) 


Integrating  Eq.  (A. 5)  with  respect  to  t '  gives  the  helical  invariant 

p  * 

P^+  r'sinte'-kQz'J-const. ,  (A. 6) 

where  p^-^mdz'/dt' ,  and  P^«rf (p^-eB^r */2c)  is  the  canonical  angular 

momentum  associated  with  the  applied  axial  field  Bq.  Subtracting 

2 

YbmVb^to“COnst*  from  Ecl*  (A.  6),  where  Ybmc  is  the  characteristic 
directed  beam  energy  and  is  the  characteristic  mean  axial  velocity 
of  the  beam  electrons,  gives 


Cb-P;+  r-  (Pi-Ybmvb)+  ~  r 'sinCe' -k0z'  )-const. , 


(A.  7) 


which  is  the  form  of  the  helical  invariant  used  in  Secs.  II  and  III 

[Eq.  (5)). 


To  determine  the  axial  invariant  associated  with  Eqs.  (A.l)- 


(A. 3) ,  we  consider  the  quantity  I.  defined  by 

z 


2  eB 


(A.  8) 


"  c!^  pefiBsin*e,“koz'*  » 

where  p'-Y'mds’/dt' ,  p'-Y’fcdr'/dt'  and  p'^Y'mt'  d0'/dt'. 
z  r  v 

Differentiating  Eq.  (A. 8)  with  respect  to  t'  and  making  use 
of  Eq.  (A. 3)  gives 
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-V‘-  & 

+c«,(6'-k0*’)r’  IT")  +  ^  “(sT 

-  8ln(e'-k0*')-p;8in(e'-k0z’)  (|£r  -k0  ||J-) 


(A.9) 


-p;co8<e'-k0r’)  (j^r  -ko  |fr)  ]  . 


Substituting  Eqs,  (A.l)  and  (A. 2)  into  Eq.  (A.9)  and  combining  terms 


gives 


(A. 10) 


which  imPlie8  that  I^-const.  is  an  exact  single-Particle  invariant. 
Adding  (eB^/ck^)^  to  Eq.  (A. 8)  implies  that 


(p‘-  5)  • 


g-  «B [P;sin (6  ’ -kQi * ) -P;cos  CO ’  -kQs ' )  ] 


(••■  StJ- 


(A. 11) 


const. , 


which  is  the  form  of  the  exact  axial  invariant  used  in  Secs.  II  and 

III  [Eq.  (6)].  For  6B  ■*  0,  we  note  from  Eq.  (A. 11)  that  C(  reduces 

to  the  axial  momentum  P'. 

z 

Finally,  subtracting  I£+2(eB0/ck0)YbmVb-const. ,  from  Eq.  (A. 4) 


gives 


-  9 

Cx"P;2+P;  +  ^  ^“VV 

-  Py 6Bcos ( 0 ' ~k0* * )+  ~  P^6Bsin(e,-k0*)-const., 


(A.  12) 


which  is  the  form  of  the  exact  PerPendicular  invariant  used  in  Secs.  II 
and  III  [Eq.  (4)]. 
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